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We propose an inductive procedure to classify A^— partite entanglement under stochastic local 
operations and classical communication (SLOCC) provided such a classification is known for A'' — 1 
qubits. The method is based upon the analysis of the coefficient matrix of the state in an arbitrary 
product basis. We illustrate this approach in detail with the well-known bi- and tripartite systems, 
obtaining as a by-product a systematic criterion to establish the entanglement class of a given pure 
state without resourcing to any entanglement measure. The general case is proved by induction, 
allowing us to find an upper bound for the number of A'^-partite entanglement classes in terms of 
the number of entanglement classes for TV — 1 qubits. 
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I. INTRODUCTION 

Entanglement resides in the root of the most surpris- 
ing quantum phenomena (of. e.g. 0). Furthermore, it 
is the main resource in the usage of quantum systems 
to process information in tasks such as cryptographic 
key distribution Q , quantum computation |3 0| , quan- 
tum state teleportation 0, quantum communication 
and dense coding |^. However a comprehensive under- 
standing of entanglement is still lacking, mainly because 
it is a highly counterintuitive feature of quantum sys- 
tems (non-separability Q) and because its analysis can 
be undertaken under different, although complementary, 
standpoints [inl |. As prominent examples the subjects of 
deciding in full generality whether a given state carries 
entanglement or not and how much entanglement the sys- 
tem should be attributed to are vivid open questions (cf. 

^3 ^'^'i references therein). This state of affairs is 
critical in multipartite systems, where most applications 
find their desired utility. 

Among others, part of the efforts are being dedicated to 
classify under diversely motivated criteria the types of 
entanglement which a multipartite system can show, ft 
is in this sense desirable, independently of these criteria, 
to have classification methods valid for any number N of 
entangled systems. One of these most celebrated criteria 
to carry out such a classification was provided in ^3 ■ In 
physical terms Diir et al. defined an entanglement class 
as the set of pure states which can be interrelated through 
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stochastic local operations and classical communications 
(SLOCC hereafter) or equivalently, as those pure states 
which can carry out the same quantum-informational 
tasks with non-null possibly different probabilities. They 
also proved the mathematical counterpart of this charac- 
terization: two states and ^' of a given system belong 
to the same entanglement class if, and only if, there exist 
invertible local operators (fLO's hereafter; that is, non- 
singular matrices), which we agree on denoting as F^^^ 
such that ^ = • • • (g) fI^I (*). Mor cover, they pro- 
vided the first classification under this criterion of tri- 
partite multiqubit entanglement, giving birth to the two 
well-known genuine entanglement triqubit classes named 
as GHZ and W classes. Later on, exploiting some acci- 
dental facts in group theory, Verstraete et al. ^3 gave 
rise to the classification of 4-qubit states. 
Regretfully none of the previous works allowed one to 
succeed in obtaining a generalizable method, fn the sec- 
ond case, the exploitation of a singular fact such as the 
isomorphism SU{2) (g) SU{2) ~ 5*0(4) is clearly useless 
in a general setting; in the first case, the use of quan- 
titative entanglement measures specifically designed for 
three qubits, as the 3— tangle to discern among dif- 
ferent entanglement classes discourages one to follow up 
the same trend, since we would have to be able to build 
more generic entanglement measures, per se a formidable 
task. However, Verstraete et al. ^3 succeeded in this ap- 
proach by introducing the so-called normal forms, namely 
those pure states such that all reduced local operators 
arc proportional to the identity matrix. These authors 
also provided a systematic, mostly numerical, construc- 
tive procedure to find the ILO's bringing an arbitrary 
pure state to a normal form. Furthermore, the use of 
these normal forms allowed them to introduce entangle- 
ment measures (entanglement monotones [T^ . indeed), 
which offered the possibility to quantify the amount of 
entanglement in the original state. In this same trend, 
other alternatives can also be found in the literature 
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[itI ITsl 01 . For completeness' sake let us recall that 
classification under SLOCC is coarser than that using 
only local unitarics. that is in which every F^''^ is uni- 
tary. Nevertheless relevant results in this realm can also 
be found in the literature |2a.21. ,.22. .23j. 
Here we offer an alternative and complementary ap- 
proach to the classification under SLOCC based on an 
analysis of the singular value decomposition (SVD) of 
the coefficient matrix of the pure state in an arbitrary 
product basis. The coefficient matrix is chosen accord- 
ing to the partition 1|2...7V with the subsequent goal 
of establishing a recursive procedure allowing one to elu- 
cidate the entanglement classes under SLOCC provided 
such a classification is known with one less qubit. The 
key feature in this scheme is the structure of the right 
singular subspace. i.e. of the subspace generated by the 
right singular vectors of the coefhcient matrix, set up ac- 
cording to the entanglement classes which its generators 
belong to. As a secondary long-term goal, the approach 
seeks possible connections to the matrix product state 
(MPS) formalism (cf. e.g. and multiple references 
therein), which is becoming increasingly ubiquitous in 
different fields such as spin chains |25|. classical simula- 
tions of quantum entangled systems j2y] , density-matrix 
renormalization group techniques and sequential gen- 
eration of entangled multiqubit states . 
We have preferred the mathematical conventions. The 
canonical orthonormal basis in will be denoted by 

{ej}j=i N (correspondingly in physics the kets \j — 

1)). Normalization is not relevant in elucidating the 
entanglement class which a state belongs to. Thus 
we will deal with unnormalized vectors and non-unit- 
determinant ILO's. In the SVD of an arbitrary matrix 
(cf. appendix), V and W will denote the left and right 
unitary matrix, whereas S will stand for the diagonal 
possibly rectangular matrix with the singular values as 
entries. In the multiqubit cases, we will agree on denot- 
ing by small Greek letters (f>,ip,... vectors belonging to 
C^, whereas capital Greek letters VPj^,... will denote a 
generic entangled state in (8) . 

The paper is organized as follows. In section |n] the en- 
tanglement of two qubits is revisited with a reformula- 
tion of the Schmidt decomposition criterion in terms of 
the singular subspaces. In section IIIII the extension to 
the three-qubit case is developed in detail and the prin- 
ciples of the generalization to multipartite and arbitrary- 
dimension systems are discussed in section Hvl We close 
with some concluding remarks in section^ An appendix 
with the relevant facts about the SVD is also included. 



II. BIPARTITE ENTANGLEMENT 

A. The Schmidt decomposition criterion revisited 

The determination of entanglement of pure states of 
bipartite systems in any dimensions, in general, and in 
two dimensions (qubits), in particular, was solved long 



ago with the aid of the well-known Schmidt decomposi- 
tion |29l l30l | , by which any bipartite state can be written 
as a biorthogonal combination 



* = 



min(Ari,A''2) 



ri=l 



(2) 



(1) 



where Ai > A2 > ■ • ■ > for all n and Ni denotes the 
dimension of subsystem i. If A„ = except for only one 
index Ai 7^ 0, then the state is a product state; on the 
contrary, if A„ 7^ for two or more indices, then the state 
is an entangled state. Furthermore, A^ coincides with the 
common eigenvalues of both reduced density operators. 
Thus, to practically determine the entangled or separable 
character of a given pure state all we must do is to com- 
pute the spectrum of pi or p2 or equivalently to analyze 
the dimensionality of their ranges. This is the backbone 
in the study of 3— partite entanglement carried out in 

Followingly in order to pave the way for a generalization 
to multipartite systems, we will reformulate the Schmidt 
decomposition criterion for bipartite systems focusing 
upon the subspace generated by the singular vectors. We 
need the next 

Definition II. 1. We will denote by 2J (resp. W) 
the subspace generated by the left (resp. right) singu- 
lar vectors, i.e. 2J — spanjwi, . . . , ffc} (resp. W — 
spanjwi, . . . ,Wk}). 

We can now state the following 

Theorem II.l. Let ^' e C™ ® C" and C(^') denote 
the matrix of coefficients of ^ in an arbitrary common 
product basis. Then ^/ is a product state if and only if 
dimSU ^ I (or alternatively dimQJ = Ij- 

Proof. Let {ei}i=i_..._„i and {fj}j=i n denote bases in 

C"" and C", respectively. Then any vector * e C™ (g) C" 
can be written as 



* = ^^c,jei(^ fj, 



(2) 



where Cy are the complex coeficients of ^E", which we ar- 
range as: 



/ Cll 



Cln 



(3) 



1 ... Cjjirt y 



The matrix C(^') = C always admits a SVD, given by 
C = ySiy^^, where V and W are unitary matrices and S 
is a diagonal matrix with entries ak (the singular values, 
indeed). Thus 



min(m,n) 
fc=l 



(4) 
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Inserting into 
and {fj}j=i,2 
composition 



and identifying new bases {ei}i=i^2 
we arrive at the well-known Sehmidt de- 



min(m,n) 



k=l 



fk- 



(5) 



The number of non-null singular values coincides with the 
rank of S, which in turn coincides with the dimensions 
of 5J and 2U (cf. appendix). □ 

From the proof we can deduce a practical method to 
recognize where a bipartite system is entangled or not: 

Corollary II. 1. Let G C" ® C" denote the state of a 
bipartite quantum system and C(^) its coefficient matrix 
in an arbitrary product basis. Then 'I' is a product state 
if, and only if, r(C(5')) = 1. 



Classification of two-qubit entanglement under 
SLOCC 



We only need one further tool to find the classifica- 
tion of bipartite entanglement under SLOCC, which is 
established as follows: 

Proposition II. 1. Let e ® denote two two- 
qubit states related by SLOCC, i.e. 



(6) 



where F^^^ and F^^^ are non-singular operators upon C^. 
Then their corresponding coefficient matrices C, C in an 
arbitrary product basis are related through 



(7) 



Proof. Just substitute \1' = J2i j=i 2 '^v^i fj ® and 
identify indices. □ 

The key idea in our analysis is to recognize the effect of 
the ILO's upon the singular vectors. If Vj (resp. Wj) 
is a left (resp. right) singular vector for the matrix coef- 
ficient C, then J't^l^fz;,,-) (resp. F^'^^\wj)) is a left (resp. 
right) "singular vector" |33 for the new matrix coefficient 
C. In order to ease the notation, we will agree here- 
after on relating * and through ^ = F'^I^ F^'^^^^'), 
which allows us to drop the transpose and Hermitian 
conjugation [sll in future considerations. 
The case of two qubits is elementary, since there is no 
much space to discuss. The bases in which the coefficient 
matrix will be expressed are the canonical orthonormal 
basis {61,62} in C^. Only two options are present: ei- 
ther dimSZn = 1 or dimSU = 2. In the first case, after 
choosing F^^l such that 



fW(z;i) = -ei, 
(Ji 



ei, 



(8a) 
(8b) 



the new coefficient matrix will turn into C = (Jo), which 
corresponds to the product state ^ ~ ei ei. We will 
agree on stating that \E' belongs to the entanglement class 
denoted by 00. 

In the second case, where cri > (J2 > 0, after choosing 
Fill and fI^I such that 



fW(wi) = — ei, 

F[21(u)i) = 61, 



fW(w2) = — 62, (9a) 

172 

F[21(w2)=62, (9b) 



( Q 5 ) , which cor- 



the new coefficient matrix will be C 
responds to the entangled state ^ = ei (8) 61 -I- 62 62. 
Now we say that ^I* belongs to the class 
The reader can readily check by simple inspection how 
in the first case the canonical matrix C has rank one, 
whereas in the second it has rank 2, as expected. In 
summary, only two classes arc possible, namely 00 and 



III. TRIPARTITE ENTANGLEMENT 

The classification of tripartite pure states is performed 
along the same lines, namely choosing the ILO's F^'l so 
that the final coefficient matrix reduces to a canonical 
one. In order to find such canonical matrices, we must be 
exhaustive in the considerations of all possibilities when 
discussing about 9J and 2U. 

The analysis of tripartite entanglement can be under- 
taken upon three possible coefficient matrices, arising 
from the three different ways to group the indices, that 
is, since * = X]ii,i2,i3=i,2 Ciii2i3eii <8) 6^^ (8) Ci^, where as 
before {ek} denotes the canonical orthonormal basis in 
C^, we have 



C^i) = Ci 



|23 



C(2) = G 



2|13 



c(3) = a 



3|12 



Cm C112 C121 C122 

6211 C212 C221 6222 

6111 C112 C211 C212 

6121 C122 C221 6222 

6111 C121 C211 C221 

C112 C122 C212 C222 



(10a) 



(10b) 



(10c) 



There is no loss of generality in choosing one of them, 
since the analysis will be exhaustive. Hereafter we will 
choose C = C^^\ Notice that now the left singular vec- 
tors of C belong to whereas the right singular vec- 
tors are in (g) C^. Also, we immediately realize that 
only two possibles options arise, namely dim21J = 1 or 
dimSJU = 2, since there are at most two positive singular 
values. The recursivity appears when classifying the dif- 
ferent structures which the subspace SB can show. The 
classification of these subspaces is performed according 
to the entanglement classes which their generators be- 
long to. In order to do that we need the following result, 
which was firstly proved in the context of entanglement 
theory in |3ll |. We offer an alternative proof in order to 
illustrate our methods. 
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Proposition III.l. Any two-dimensional subspace in 
contains at least one product vector. 

Proof. Let be a two-dimensional subspace of ® . 
With no loss of generality two entangled vectors can be 
chosen as generators of V with coefficient matrices given 
by Ci = I and C2 being an arbitrary rank-2 matrix in the 
product canonical basis. Then it is always possible to find 
non-null complex numbers a and f3 such that al + PC2 
has rank one|39j. □ 

In other words, this proposition shows that 
spanl^*!, ^'2} always equals either span{(/)i 02 '8''!/'2} 
or span{0 ® V'j^}: where implicit are the assumptions 
that different indices denote linear independence and 
in the last case only one product unit vector can be 
found. Thus, with the same convention, the right 
singular subspace 211 can show six different structures, 
namely span{(^ ® -0}, span{4'}, span{0 ® i'i:<j) <8> ^"2}, 
span{0i V:02 ^ V'}: span{0i -01 , 02 ® V'2} and 
span{0 ® 0!, VE*}. We pursue ProD. HTLll a step further: 

Proposition III. 2. Let W he a two-dimensional sub- 
space in ® C^. Then 20 = span{0Cx) 1^, xp} if and only 



if, 2IJ = span{0(8"^, 00</3-|-00</3}, where denotes linear 
independence. 

Proof. Suppose (g) is the only product vector in 21J 
(up to normalization factors) . Its orthogonal vector in W 
will be an entangled vector with coordinates in a product 
basis {0,0}(8){v?, (f} given by (0, ai2, a2i, "22), i-c it will 
be of the form 0(g)(p-|-0®v+a0®i^, with a £ C. Since 0® 
if must be the only product vector in SB, it necessarily has 
to be a = 0; otherwise could it always be possible to find 
a, /3 £ C such that q0 (S^ip-\-(3(^(t)®(p-\-(j>^(p-\-a(t>®(p) 
is another product vector (/? = aa) . 
Suppose now that 211= {0®(/?, 0(g)(^-|-0®iy9}, then 
a0 ®(y3-|-/3(0(g)(^-|-0®iy9) is a product vector if, and 
only if, /3 = 0, i.e. if it is the original cj) ^ (p. □ 

We can now state our result, already contained in 
with different criteria: 

Theorem III.l. Let * e ® ® be the pure state 
of a tripartite system. Then 5" can be reduced through 
SLOCC to one of the following six states, which corre- 
sponds to the six possible entanglement classes, according 
to the following table: 



Class 


Canonical vector 


Canonical matrix 


w 


000 


ei (8) ei (g) ei 




^10 0^ 

^0 




span{0 g) -0} 




ei (g) ei ei -1- ei (g) 62 (g) 62 




^ 1^ 

^0 oy 




span{vl/} 




ei (g) ei ei + 62 g) ei (g) 62 




^10 0^ 

^0 1 oy 




0(gC2 


03^-^2 


ei (g ei g) ei -f 62 (g) 62 gi 61 




^10 0^ 

^0 1 oy 




g)0; 


GHZ 


61 g) 61 g) 61 -1- 62 g) 62 g) 62 




^10 0^ 
^0 ly 




span{0i g) -01, 02 g) ■02} 


W 


61 g) 61 g) 62 -f 61 g) 62 g) 61 -)- 62 gi 61 g) 61 




^0 110^ 

^1 oy 




span{0i -01, ^E'} 



Proof. We discuss depending on 211: 



1. 2IJ = span{0g)-0}. In this case, wi 
the ILO's , fc = 1, 2, 3 so that 



1 

0-1 

ei, 

61. 



-ei, 



igj^j. Choose 



(11a) 



(lib) 
(11c) 



Then the new coefficient matrix will be 

/l o\ 



C 

















■0 











7 



1000' 
,0000, 



(12) 



which corresponds to the state 61 g) 61 g) 61, and 
where the dots • indicates the irrelevant character 
of that entry. 



5 



2. 2IT = span{^'}. In this case wi = ® t/ji + 02 (8) '02 ■ 
Choose the ILO's so that 



wi and W2 are hnear independent (orthonormal, 
indeed). Choose the ILO's so that 



0-1 

F^^H^Jl)=ei, F[31(V;2) = 62. 

Then the new coefficient matrix wih be 



(13a) 

(13b) 
(13c) 



C 



(Ti 
, 0, 



/l l\ 

[■■■■) 



noon 
yo oy ' 

which corresponds to the state ei 

62 (E) 62. 



(14) 



6i (X" 6i + 6i 



3. 2U = (8) C = span{0 (g) 0i, '02}- In this case 
wi = ^110 Vi + Mi20 ^ V'2 and W2 = Ai2i0 ® "01 + 
/i220®V'2, where the matrix [/ly ] has rank 2. since 
wi and i(;2 are linear independent (orthonormal, 
indeed). Choose the ILO's so that 



i^i'' = [^^r(ei)^^r(e2)] = K-]-\ 
F[2l(0) = ei, 

F[3l(V.l) =61, F[31(V^2) =62. 

Then the new coefficient matrix will be 



(15a) 

(15b) 

(15c) 
(15d) 
(15e) 



C 



Mil Ml2 
vM21 M22, 



^ 

1 




/^Mn Mt2 0\ 

M21 M22 



M 0\ 

yo 1 oy ' 

which corresponds to the state 61 
61 ® 62. 



(16) 



61 (XI 61 



62 ' 



4. W = C"^ (E) ip ^ span{0i O -0? </'2 ® "0}- In this case 

Wi = ^1101 V + M1202 f?) -0 and W2 = M2101 ® + 

/i2202 ® "07 where the matrix [^y ] has rank 2, since 



f]^1(i;i) 



-61, 



Ff'(^^2) 



762, 



^^r = N''(ei)F^^'(62)] = K-] 



62, 



F[2l(0i) = ei, F[21(02) 
F[31(V') = 61. 
Then the new coefficient matrix will be 



(17a) 

(17b) 

(17c) 
(17d) 
(17c) 



C = 



Mil M12 
vM2i M22 , 





M21 M22 



V 



7 



1000 
,0010, 



(18) 



which corresponds to the state 61 
62 (8) 61. 



61 (Ki 61 



■ 62 ' 



5. S2U = span{0i -01,02 (8) V'2}- In this case Wi = 
^1101 (8) -01 + Aii202 8) ^"2 and W2 = /^2i0i 8) V'l + 
Ai22 02 ® "02, where the matrix [/^i^] has rank 2, since 
wi and W2 are linear independent (orthonormal, 
indeed). Choose the ILO's so that 



FI'\v2) 



FW = [Fi^'(6i)Fi^'(62)] = [/iy 
FW =FWFfl, 
F[2l(0i) = ei, F[21(02) = 62, 
F[31(V;i) = 61, ^[■■'1(V'2) =62. 

Then the new coefficient matrix will be 



C2 
1-1 



62, 



(19a) 

(19b) 

(19c) 
(19d) 
(19e) 



C = 



Mil M12 
vM2i M22, 





/mu ^il2\ 
M21 ^^1 



V .... y 

/i o\ 
1^0 ly ' 

which corresponds to the state 61 (8 61 (8> ei + 62 

62 62. 



(20) 
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2IT = span{t/)i (Eiipi, In this notation, remember 
that imphcit is the assumption that only one prod- 
uct unit vector can be found in 2n. In this case 
'I' can be chosen so that 5" = ® ■(/'2 + 02 ® V'l 

Thus the 



(this is the statement in Prop. IIII.2|) 
singular vectors can always be expressed as wi = 
Mil {4>i ® "02 + 02 V'l ) + Mi2</'i <8> ipi and W2 = 
M21 {4>i €5 ■02 + 02 <8) V'l ) + M2201 ® '01 J where the 
matrix [jJ^ij] has rank 2, since wi and W2 are lin- 
ear independent (orthonormal, indeed) . Choose the 
ILO's so that 



-ei, 



re2, 



FW = [FW(ei)FW(e2)]-Kr\ 

F[2l(0i) = ei F[2l(02) = e2, 
i^[3](^i) = ei F[3l(^2) = e2. 

Then the new coefficient matrix will be 



(21a) 

(21b) 

(21c) 
(21d) 
(21e) 



C 




^/^21 M22 / 

fai 0^ 
\ (72 0^ 

^0 1 1 o\ 
,10 0/' 



which corresponds to the state ei 

62 (8 Ci -f 62 ® ei (g) Ci. 



/mI2 /^ii Mn 0\ 

M22 ^^21 M2I " 



V 



7 



(22) 



'62 + 61 



Since there is no more options for the subspacc S2IJ we 
have already considered all possible alternatives. 

□ 



In conclusion, we have found that there are six classes 
of entanglement, named after as 000, Oi^^'^j^, GHZ 
and W . The theorem also indicates how to practically 
classify a given state compute the SVD of its coeffi- 
cient matrix and elucidate the structure of spanjwi, W2}- 
We include a further proposition comprising the practi- 
cal implementation of this result. We need to introduce 
the following definition. 

Definition III.l. Letwj = ei^Wji+e2<E)'Wj2 € C^(g)C^ 
be an arbitrary vector. We associate a two-dimensional 
matrix Wj to Wj by defining 



Uji Wj2\. 



(23) 



This definition will be mainly applied to the right sin- 
gular vectors of the coefficient matrix C. As usual, the 



singular values of C will be denoted by ak , in nonincreas- 
ing order, and (j(A) denotes the spectrum of a matrix A. 
Our proposal to implement the preceding result is 

Tiieorem III. 2. Let "if denote the pure state of a tri- 
partite system and C*-*' its coefficient matrix according 
to the partitions i\jk (cf. H10a|) - (|10cf) ). Then 

1. 'J belongs to the 000 class if, and only if, v{C^^^) = 
1 for all i = 1,2,3. 

2. 5* belongs to the ^1^23 class if, and only if, 
r(C(i') 1 andr(C('^)) ^ 2 for k =^ 2,3. 

3. ^ belongs to the 02*1' 1*3 class if, and only if, 
r(C(2)) = 1 andr(C('=)) ^ 2 for k = 1,3. 

4- belongs to the 03'^^2 class if, and only if, 
r(C(3)) = 1 and r(CW) ^2 fork 2. 

5. 'J belongs to the GHZ class if, and only if, one of 
the following situations occurs: 

i. r(CW) = 2 for all i = 1,2,3 and t{Wi) = 
r(W2) = 1. 

a. r(CW) = 2 for all i = 1,2,3, r{Wi) = 2, 
r{W2) = 1 and a{W{^W2) is non-degenerate. 

ill. r(CW) = 2 for all i = 1,2,3, r(M^2) = 2, 
v{Wi) = 1 and g{W^'^Wi) is non-degenerate. 

iv. r(CW) = 2 for all i = 1,2,3, t(Wi) = 2, 
r{W2) = 2 and a{W^^W2) is non-degenerate. 

6. belongs to the W class if, and only if, one of the 
following situations occurs: 

I. r(C«) = 2 for all i = 1,2,3, y(Wi) = 2, 
r(W^2) = 1 and cr{Wi^W2) is degenerate. 

a. r(CW) = 2 for all i = 1,2,3, t(Wi) = 1, 
^{W2) ~ 2 and a(W2^Wi) is degenerate. 

ill. r(C«) = 2 for all i = 1,2,3, r{Wi) = 2, 
r(W2) = 2 and a{Wi^W2) is degenerate. 

Proof. We will exclusively concentrate upon the suf- 
ficiency, since the necessity directly follows from the 
canonical form of each class. 

The first four cases are elementary, since it is a matter 
of detection of the vector which factorizes. The final two 
cases correspond to true tripartite entangled states. If 
^{Wk) = 1 ior k — 1,2, it is clear that there exist two 
product vectors belonging to W, thus 5* belongs to the 
GHZ class. If t{Wi) = 2 and r{W2) = 1 we need to check 
whether an ILO applied upon the first qubit can reduce 
the rank of the transformed Wi. As it can be deduced 
from the preceding proofs, an ILO upon the first qubit 
amounts to constructing a linear combination between 
the two right singular vectors, which is equivalent to find 



new matrices Wj = F^^Wi 



F^fW2, with j 



1,2. If 
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r(Tyi) = 2, then by multiplying this expression to the left 
by W^^, we have 

FWl + F^T^fVa. (24) 

It is immediate to realize that it is possible to reduce 

the rank of Wi to 1 and to choose f}]'^ such that F^'^^ 

''J 

is nonsingular provided the spectrum of W^'^W2 IS non- 
degenerate, in which case \1/ belongs to the GHZ class. 
If the spectrum is degenerate, thus both eigenvalues be- 
ing null, no further reduction is possible and ^E* belongs 
to the W class. The symmetric case runs along parallel 
lines. 

Finally if r(iyi) = r(W2) = 2, reasoning along similar 
lines if both eigenvalues of W^^W2 are equal, only one 
rank can be reduced keeping the nonsingularity of i^I^l 
and 'i> belongs again to the W class, whereas if the eigen- 
values are different, both ranks can be reduced to 1 keep- 
ing the nonsingularity of F^^l and belongs to the GHZ 
class. □ 

As a final remark let us indicate how close, despite 
the apparent differences in the a ppr oach, our analysis 
runs parallel to that performed in [IjI : the ranges of the 
reduced density operators are indeed generated by the 
corresponding singular vectors, and the study of these 
ranges drove them and has driven us to the same final 
result. The change of method is motivated by the at- 
tempt to find a generalizable criterion not using entan- 
glement measures specifically built upon the number of 
qubits of the system, such as the 3— tangle With 
this approach it is not necessary to consider at any stage 
the reduced density matrices and entanglement measures 
upon them. A strongly related approach can be found in 
[3II33. 

IV. GENERALIZATIONS (A > 4) 

The generalization of the preceding approach to pure 
states of arbitrary multipartite systems is two-folded. On 
one hand, the generalization to multiqubit states can be 
implemented inductively: 

Theorem IV. 1. // the entanglement classes under 
SLOCC are known for N qubits, the corresponding en- 
tanglement classes for N + \ qubits are also known. 

Proof. Wc proceed by induction. We have proved in pre- 
ceding sections that this statement is true for N = 2 and 
have explicitly found the entanglement classes for A^ = 3. 



For a given {N + l)-qubit system, write the coefficient 
matrix Ci|2...Ar+i = C. Because of the induction hypoth- 
esis one knows in advance the classification of the right 
singular subspaces of C according to 211 = spanj^'i} if 
dim an = 1 and 2U = spaufvlf^, vj/^ } if dim 22} = 2, where 
each and belong to one (possibly the same) of 
the entanglement classes of A^ qubits. Choose the ILO's 
i?[2] (g) . . . (g) so that the two first columns of W 

(the transformed right singular vectors) are expressed as 
linear combinations of the canonical vectors of the en- 
tanglement classes corresponding to the structure of 20 
and choose the ILO F^^^ so that ySM^^ drops out as 
many non-null entries as possible (typically F^^^ will be 
the inverse of a rank— 2 submatrix of VF^). The result is 
the canonical matrix for an entanglement class of A^ + 1 
qubits. □ 

There is an important remark in the preceding induc- 
tive construction, already stated in |l2l | and explicitly 
shown in |l3l| : there will be a continuous range of states 
with a similar right singular subspace but with no ILO's 
connecting them. Let us illustrate this peculiar fact with 
an explicit example. When considering 4— partite entan- 
glement, there will exist 45 a priori structures of the 
right singular subspace of the coefficient matrix, aris- 
ing from 6 possible one-dimensional right singular sub- 
spaces 2n = span{^'}, where belongs to one of the six 
entanglement classes of A^ = 3, times 4 possible sites 
for the fourth added qubit, plus 21 possible bidimen- 
sional right singular subspaces 20 = spanj^I'i, vj/j}, cor- 
responding to the (^^2~^) ways to choose the classes for 
A^ = 3 which and '^2 belong to. An example will be 
2IT = spanjOOO, GHZ}, with the already convention that 
only one product vector and no Oi'^jk belongs to 211, i.e. 
2U = span{0i ig) ^pi ® ipi,4>2 'S> ip2 'S' 4'2 + 4>2 <f2 V'2}, 
where the vectors with ~ are pairwise linearly indepen- 
dent. In order to only have one product vector and the 
rest being GHZ vectors, we must have 113 (up to per- 
mutations) 22J = span{0(8)(^(g) V'', 0i8)</'0V' + <^<55<^(8 V;}, 
with ip' 7^ tp, ijj. 
Recalling that 

Wi = ® (p ® Ip' + ^112 {4> ® (fi ® ip + (j) ® (p ® ^) , 

(25a) 

W2 = /i2i0 <8) (8) V'' + M22 ® V? ® "0 + ® ® Vj) , 

(25b) 

where the matrix [/iy] = (^Jj^ ) will be non-singular, 
it is immediate to find ILO's i^[2] ^ i^[3] ^ g^^^ that 



g, ^[3] ^ F^'^\wi) = finei 62 ® V + M12 (ei «) ei «) d 62 ® 62 «> 62) , 

g, ^[3] ^ ^[4] (y^^-) ^ ^^^g^ 52 g) + ^22 (ei (g) ei ® ei 62 62 «) 62) 



(26a) 
(26b) 
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which corresponds to a coefficient matrix given by 



/mi2 /inVi MiiV'2 /Xi2\ 

M22 pi2lV'l ^J*2l^2 /X^2 

C = VY1 , (27) 



V 



5x8 



of the transformed in the canonical basis. Choosing 
Fill so that 



we arrive at 



C 



VI V'2 

,10 00001, 



where the coefficients ij^i corresponds to the coordinates which corresponds to the canonical vector 
I 



(28) 



(29) 



ei (8) ei €5 62 (8) ■0* + 62 ® ei (g) ei ei + 62 62 ® 62 ® 62 = (-0* 7^ ei, 62) 

= 10010*) + 11000) + lllll) (IV/) ^ |0), |1)) 



(30) 



Thus, different ^ will yield different entanglement classes 
under non-singular local operators F^^l ® ■■■ ® F^^\ 
Notice that this vector belongs neither to the GHZ4 
class nor to the W4 class nor to the $4 class (containing 
the cluster state of four qubits -see below). It is a 
peculiar feature that two infinitesimally close states 
could belong to distinct entanglement classes, so a 
deeper elucidation of this point is on due and will be 
carried out also elsewhere. For the time being, we will 
agree on attributing all states reducible to H30() by ILO's 
fI^I ® ■ • ■ ® fW, independently of the particular vector 
Tp, the same entanglement properties under SLOCC and 
analogously for arbitrary iV-partite multiqubit systems. 
This allows us to find an upper bound for the num- 
ber of genuine (TV -|- l)-partite entanglement classes. 
Firstly, notice that e.g. the right singular subspace 
20 = spanjOOO, 000} in the 4— partite case actu- 
ally contains structures with different properties, 
namely [ilj 22J = (8 (and permutations), 
2U = span{0 8) Vi > <8"^2 <8i V'2} (and permutations) 
and 22J = span{0i ^ ipi ^ 0i,02 (8 <P2 0$ ^'2}- AH of 
them drives us to at least one factor qubit in the final 
canonical state, except one, that is, there will correspond 
one right singular subspace structure spanj^fi, ^'2} to 
each genuine {N + 1)— entanglement class. 

This is rigorously proved in the following 

Proposition IV. 1. Let Wn be the right singular sub- 
space of the coefficient matrix in an arbitrary product 
basis of an N -qubit pure state. If TOn is supported in 
a product space Wm = ^ ® then the state be- 

longs to a product class , where denotes a class of 
(N — I) -partite entanglement. 

Proof. Under the above assumption, Wj = ip ® Wj, j = 
1,2, with G and wj e C^^^"^^ We can always find 
an ILO fI^I such that 

Wj ^ 61 (g) Wj, (31) 



where also Wj g C^^^"^^ hence Wn = En ® Wn-i, 
where En denotes the Weyl matrix En = |6i)(6i|. Since 
we can always write Eat = Fh ® S^r-i, the coefficient 
matrix can always be written as 

Cn = VJ:nW^ ^V{En<E)J:N-i){En<E)WN^iy 
= En ® {v^N-iWl_^^ . (32) 

The remaining ILO's F^^l and F^^\ j > 2, can always be 
chosen so that 

CN = En®CN-i, (33) 

where Cn-i denotes a canonical matrix of an [N — 1)- 
partite entanglement class. This proves that the second 
qubit factorizes, as the reader may check. □ 

With appropiate permutations, this result applies to 
any qubit. If we denote by M{N) the number of 
A^— partite entanglement classes, there will be at most 

^ + 11 (34) 

genuine entanglement classes for + 1 qubits. Besides, 
the number of degenerate (N -\- 1)— entanglement classes 
will be at most {N -I- 1) x M{N) (corresponding to the 
1 possible factor positions which the (A^-|- l)th qubit 
can occupy), thus 

Corollary IV. 1. Let M{N) denote the number of 
N— partite entanglement classes under SLOCC. Then 

M{N + 1) < ^M(iV) [M{N) + 2A^ + 3] . (35) 

The equality will be in general unattainable, since, as 
in the case of tripartite entanglement, only a few distinct 
true entanglement classes exist, coming out from the only 
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actually different structures which the right singular sub- 
space can adopt (only two in the case of tripartite sys- 
tems; cf. proposition IIII. 1 1 . 

Let us call reader's attention on the fact that these re- 
sults allow us to view all state space of N qubits divided 
into blocks, each one parametrized by a right singular 
subspace structure and corresponding to our broad-scnsc 
entanglement classes, and within which the difference be- 
tween states stems from a (possibly several) continuous 
parameter. The number of these blocks for N + 1 qubits 
is upperly bounded by the recursive relation H35|) . 
Another benefit of the present approach arises when de- 
ciding whether two states belong to the same entangle- 
ment class or not. This is stated as a corollary: 

Corollary IV.2. Let $ G (C^)*^^. Let and 211$ 
be their respective right singular subspaces. Then a nec- 
essary and sufficient condition for $ to belong to the 
same broad-sense entanglement class under SLOCC is 
that 2U* and 2U$ have the same structure, i.e. that they 
are generated by entanglement- equivalent vectors. 

Proof. The result follows immediately both from con- 
struction and from the convention on the definition of 
the broad-sense entanglement classes. □ 

As an example, let us include a one-line 
proof that the 4-qubit GHZ state \GHZi) = 
^ (|0000) -I- |1111)) and the cluster state |04) = 

i (10000) + 10011) Hj_[1100) - |1111)) do not belong 
to the same class [83 ■ Their respective right singular 
subspaces are WcHZi = spanjei ei (g) ei, 62 62 ® 62} 
and SIJ^j = spanjei ® 5'"'', 62 ® ^~}, where denote 
two-qubit Bell states. It is immediate to conclude that 
they are different, since none Cj ® Cj ® Cj belong to 
2f 04 (write the coefficient matrix of a generic vector in 
in terms of two coordinates a and /3 and check 
that it is impossible to choose the latter so that the 
matrix corresponds to ej (g) g ej). These states belong 
to the respective so-called GHZ4 and $4 classes, 
characterized by the above right singular subspaces. 
On the other hand, to find a wider generalization one 
can focus upon arbitrary dimensional entangled systems. 
The leit motiv is still the same, with the important 
exception that the dimension of the right singular 
subspace can grow up to the dimension of the Hilbert 
space of the first subsystem. Thus the analysis of the 
possible structures which SB may adopt is now much 
more complex. 

We include as an illustrative immediate example the 
analysis of all entanglement classes under SLOCC of 
any bipartite (iVi x Af2)-dimensional system: there exist 
min(A^i, A^2) entanglement classes, which can be denoted 
as 00 = ^-j^, ^-J, . . . , 'iftun{N,,N2)' '^^o^^ canonical 
states will elementarily be ^ ^i' ^o^' each class 

They correspond to canonical matrices given by 
J2i=i \^\^^^){^\^^\^ so that we can state the following 



Theorem IV.2. Let * G C^i (giC^^ be the pure state of 
a bipartite quantum system with coefficient matrix in an 
arbitrary product basis denoted by C(^). Then ^' belongs 
to the class, k = 1, 2, . . . , min(A^i, Af2)7 */; o-nd only 
if r(C(1')) = fc = dim5J = dim2IJ. 

Proof Let QJ = span{^fc}fc=i^... .„<n,in(Ar,,Ar2) and 2IJ = 

span{(^A;}fc=i,... ,n<min(Afi,Ar2)- Choosc F^^^ and fI^I so 
that 

F^'\M = —ek, (36) 

F^'Kvk) = efc. (37) 

Then the coefficient matrix (in blocks) will turn out to 
be 

C^i ^" '^^^-'^ I . (38) 

□ 

For more general cases, the difference stems solely in 
the higher computational complexity. 

V. CONCLUSIONS 

We have developed a recursive inductive criterion to 
classify entanglement under SLOCC in multipartite sys- 
tems in pure states which allows one to find the entangle- 
ment classes for -I- 1 qubits provided this classification 
is known for A^ qubits. The method rests on the analy- 
sis of the right singular subspace of their coefficient ma- 
trix, which is chosen according to the partition 1|2 . . . A^, 
hence a 2 x 2-^^"^ rectangular matrix. Then one must elu- 
cidate the classification of the one- and two-dimensional 
right singular subspaces according to the entanglement 
classes which their generators belong to. As a conse- 
quence, this construction reveals a systematic way to 
detect the entanglement class of a given state without 
resorting to quantitative measures of entanglement. In 
arbitrary-dimensional generalizations, the same scheme 
must be followed with the exception that the dimension 
of the right singular subspaces is higher and their struc- 
ture now depends on several generators. 
For A^ > 4 it has been showed that within each right sin- 
gular subspace structure, there could exist a continuous 
infinity of states not connected through invertible local 
operators. Additionally, up to this continuous degree of 
freedom within each right singular subspace structure, 
we have found an upper bound for the number of classes 
on A^ + 1 qubits in terms of the number of classes of A^ 
qubits. 

As a final remark, let us conjecture that a possible con- 
nection with the MPS formalism is probable to exist. In 
this formalism (cf. and multiple references therein) 
any pure state is written as 

ii ...ijv 
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so that adjoining a further {N + l)-th qubit amounts 

to adjoining a further aJ^'^Y' rnatrix in the trace giving 
the coefficients. In the analysis carried out above, this 
last added qubit is equivalent to increase the dimension 
of the right singular subspace dimW^ — * dimSUAr+i = 
2 X dimWiy. Our conjecture is that the structure of Wjy 
should be read from the properties of the N matrices 
A^^''\ so that the succesion of structures of SUat should 
run parallel to that of the matrices A^^^^ , . . . , A^"' . 

APPENDIX A: THE SINGULAR VALUE 
DECOMPOSITION 

We include the relevant properties of the SVD of an 
arbitrary matrix and suggest the interested reader to 
consult e.g. [s^ for a comprehensive analysis of this de- 
composition with the corresponding proofs. The set of 
m X n complex matrices will be denoted as usual by 
■M.m,n(jC) = Aim.n and the group of unitary matrices 
of dimension k will be denoted by U (fc). The main result 
can be stated as 

Theorem A.l. (Singular Value Decomposition) 

Let Q G A^m,n- Then Q can always be decomposed as 



singular vectors of Q. The positive entries of Yi are the 
singular values of Q. 

Notice that with this definition any mx n dimensional 
matrix will have m left singular vectors and n right sin- 
gular vectors; since the relevant singular vectors will be 
those associated to non-null singular values, we agree, as 
usual, on referring as singular vectors only to the latter, 
i.e. to those Vk and Wk for which ak > 0. Another com- 
mon convention is the decreasing order of the singular 
values in the diagonal of E: ai > a-2 > ■ ■ ■ > 0. 
The singular vectors are highly nonuniquc or cquivalently 
there always exist another unitary matrices V and W 
such that Q = VT,W^ , where these new unitary matrices 
depend of the former V and W and the multiplicities of 
each singular value [s^. However this fact has not been 
exploited in the text. 

One of the main consequences of the SVD is that the 
rank of a given matrix Q coincides with the rank of E, 
i.e. with the number of positive singular values, which, in 
turn, coincides with the dimension of the subspace gen- 
erated by the left (or right) singular vectors. This is the 
basis to the analysis of entanglement of a pure state upon 
its coefficient matrix in a product basis performed in the 
text. 



g = VT.W\ (Al) 

where V G U{m), W G U{n) and E G M^.n is a diagonal 
matrix with non-negative entries, i.e. Ynj = cridij, with 
i = 1, . . . ,m, j — 1, . . . ,n and ak > for all k. 

The columns of V and W and the positive entries of E 
receive a special name: 

Definition A.l. The columns of V ^ [vi V2 ... Vm] 

(resp. W = [wi W2 ■-. Wn]) are the left (resp. right) 
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